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Quarterly Progress Report, Contract No. NAS 5- 27 483 

Landsat-D Investigations in Show Hydrology 
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terrain. Tenth International Symposium, Machine Processing of Remotely Sensed 
Data, 230-234. 

Dozier. J., 1984, Snow reflectance from Landsat-4 Thematic Mapper, IEEE TYans actions 
on Geoscience and Remote Sensing. G&22, 323-328. 

Frew, J., 1984, Registering Thematic Mapper imagery to digital elevation models. Tenth 
International Symposium, Machine Processing of Remotely Sensed Data, 432-435. 

(Copies attached.) 

* Development at Twes tr aam Model (paper in preparation) 

INTRODUCTION 

Two-stream methods provide rapid approximate calculations of radiative transfer in 
scattering and absorbing media. Although they provide information on fluxes only, and not on 
intensities, their speed makes them attractive alternatives to more precise methods. 

Meador and Weaver [1960] provide a comprehensive, unified review of two-stream methods 
for a homogeneous layer, and solve the equations for reflectance and transmittance for a homo- 
geneous layer over a non-reflecting surface. In this paper, we show how any of the basic kernels 
they describe for a single layer can he extended to a vertically inhomogeneous medium over a 
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surface whose reflectance properties vary with illumination angle, as long as the medium can be 
subdivided into homogeneous layers. The method never produces nonphysical results (such as 
negative albedo), and it is computationally stable, without any restrictions on the optical proper- 
ties of the layer. . 

The outline of the procedure is: 

1. For the bottom layer, over a surface whose reflectance properties are known, solve the 
two-stream equations to calculate reflectance and transmittance for direct and diffuse irra- 
diance. 

2. The bottom layer can now be considered the new "surface," whose reflectance properties 
are known, so successively calculate reflectances and tranamittances for each layer to the 
top of the medium. 

3. Now start at the top with known incident flux, and successively calculate upward and down- 
ward fluxes for each layer. 

TVO-5TREAM APPROXniATIQN FOR A BOHOGBfEOUS LAYER 

Uaador and Wtavtr [1980] review the derivation for two-stream methods, so the treatment 
here is brief. They show that all of the common variants, including the Eddington and delta- 
Eddington methods, are different cases of the same differential equations. 

The optical properties of a layer are described by: 

0 Single-scattering albedo, the portion of the radiation incident on a scattering element 
that is reflected, instead of absorbed. It is defined by the ratio of scattering to extinc- 
tion efficiency. 

g Mean value of the cosine of the scattering angle, measured from the forward direction. 
g= + l corresponds to completely forward scattering, -1 to completely backward 
scattering, and 0 to isotropic scattering. 

t 0 Optical thickness of the layer. 

Direct flux Eq, measured perpendicular to the beam, is incident at angle arccos/xo- The 
diffuse upward and downward fluxes F* and F~ at optical depth t within the layer are described 
by the following pair of equations for any two-stream approximation: 

7i** - 7 *r- -*Eot*^ /Mo (i) 

^~= 7.^* - 7t r- + o»ffo 74 t" T ' ,<D (2) 

The layer is homogeneous if the y's are independent of t. The choice of the 7*3 is deter- 
mined by the particular approximation to the scattering phase function and the radiation inten- 
sity distribution. One constraint is that 7 s + 74 = 1 , because of energy conservation. Idtador and 
Wtavar [1980] give expressions for the 7 's for 7 different two-stream approximations. 

GENERAL SOLUTIONS 

General solutions to ( 1 , 2) differ for non-conservative (o»< 1 ) and conservative ( 0 = 1 ) 
scattering [fiscombt, 1977]. 


1. Nan-Conservative Scattering 

A general solution for the case 0 < 1 , is 


= E 0 
F~ =* E 0 


Ci(t + 7i) t ,r C«(l-7 i) t -fr _ 


7« 


7« 


Ci eh 1 +• C,e-* r - P a*** 


(3) 

(4) 


- * ** 


ORIGIrtAl PagE 13 
OF POOR QUALITY 


where 

f 

P-kn o(ai/4o+7«) 

9 “ *A<o(«^o-7t) 

ai = 7i74+7i7S y 

a* = 7 i7s+7i74 

* = S 

1-<W 

Ci and C, are arbitrary constants. Given top and bottom boundary conditions, their expressions 
can be derived, and the general solution can be used to find reflectance and transmittance of a 
layer. 


2. Co nser vative Scattering 

When w~ 1, a general solution is 

F* = £„[c, + C,t— / i 0 (7iMo-7a)«" T/Mo j 


Ci + Cb(t- - Po(7iAc + 7*) • _T/ *‘ q 

7i 


r = E 0 

where Ci and C| again are arbitrary constants. 


(5) 

( 6 ) 


BOUNDARY CONDITIONS 

The coefficients Ci,C B depend on the boundary conditions. The usual top boundary condi- 
tion is that there is no diffuse irradiance at the top of the layer. 

r-(o) = 0 (7) 

In atmospheric radiation models the usual bottom boundary condition is simple. Unador 
and Weaver [1960] for example use a flat, unobstructed, zero-albedo surface; Wiscomb* and War- 
ren [1980] use a flat, unobstructed, Lambertian surface. To apply an atmospheric model over 
the land however, we must allow for the surface reflectance properties to vary with incidence 
angle, and for irregular terrain such that the illumination angle at the surface is not necessarily 
the same as at the top of the atmosphere. Moreover we must recognize that a significant portion 
of the incident diffuse irradiance comes from reflection from adjacent terrain. 

Define the following quantities; 

/it cosine of illumination angle at the local surface; 

p'tf bihemispherical reflectance of the surface; 

p’ t (p,) directional-hemispherical reflectance of the surface at incidence angle cos~V«; 

Vg terrain “view factor” for incident diffuse irradiance; 

V, terrain view factor for incident direct irradiance. 

Then the lower boundary condition is 

^(To) * *o«" V *[p , .N.+W'e , IS] + r*(T 0 ) Ip's [1-(1 -p's) V 4 ]\ (8) 

The terms bihemispherical and directional-hemispherical reflectance are defined by 
Nicodamus et ai. [1977] to describe the hemispherically integrated reflectance to diffuse and 
direct irradiance. The view factor V 4 for diffuse irradiance accounts for the portion of the hemi- 
sphere “seen” by a point, weighted by the angles between the point and the surrounding terrain 
in the range of azimuth directions. S is the slope angle, E is the exposure, the direction that 
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the slope faces, and H f is the angle to the horizon in direction p. 

w 

V 4 - J [ cos 5 cos H f + sin S camH f ccm(y-E) ]*dp (9) 

The view factor V, for direct irradlance is similar to V 4 , except that the integrand for each direc- 
tion p is weighted by the cosine of the illumination angle on the slope to the horizon, it. by 
[>uocos// f - (l-zio^sin//, cos(pg-p)], where this te.-m is set to zero if it is negative, po is the 
azimuth of illumination. 

Simplified bottom boundary conditions can be implemented by modifying some of the above 
terms. If the surface is Lambertian, for example, set p' t =pV If the surface is flat, set /i, = Po- 
If the surface is unobstructed by surrounding terrain, set V 4 = V, = 0. Even in rugged terrain we 
can often assume V, V 4 . 


museums or a homogeneous layer 


Directional-hemispherical reflectance of the layer (also called direct albedo) is the upward 
radiance emerging from the layer, integrated over the hemisphere, divided by the direct 
incident irradiance. Bihemispherical reflectance (also called diffuse albedo) is the upward radi- 
ance divided by the downward diffuse irradiance, and can be derived by integrating directional- 
hemispherical reflectance over the hemisphere. 

p. Gio) =- (10) 


PS 


= £10 1 
^(0) 



(11) 


Direct albedo is calculated by finding the values for C*i and C a in the general solution, given 
the boundary conditions, and then substituting these values in (10). 

Integration of p v over /a. to obtain p a is done numerically, although for particular approxi- 
mations to the phase function analytic integration is passible. For example, Wiscomb* and War- 
ran [1980] give an analytic integration for the delta-Eddington approximation. 

1 . Non -Consul native Scattering 

For u< 1, direct reflectance is 

P.M - 

♦ (12) 


where 

D* = tTs" 1 *® + G* # |T ° 

= 7*±7i 

6 *P*s[l-(l-p’s)Ks] 

if= -2 — As[d(o^o+7 4 )-a«Po + 7s] 
Pc 

7) =P'iP* +P'S v. 


When the layer is semi-infinite (i.e. when To 

_ (v- >/. x _ ® (7i-t74) 

"• <Po) ‘7.<i+<*0 


®) reflectance is: 


( 13 ) 
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2. Cun — i » aU — Scattering 

In the conservative scatttering caw, when a = 1, direct reflectance la 
m l - ±-x 

|(l-d) (jtjjQ+ji) ~ e“ V> * 0 1(1 -d) ( 71 P 0 + 7 O + 7 }— 1| 1 


where 

D, = l + ro(7,(l-d) 



TBAN3HTTANCE OF A HOMOGENEOUS LAYER 

We distinguish three kinds of transmittances: 

Directional transmittance T, is the proportion of the direct irradiance incident on the top 
of the layer that is transmitted without scattering or absorptir.n. 

r.U) = 0~ r ° /M ° (15) 


Directional-hemispherical transmittance T* is the proportion of the direot irradiance 
incident on the top of the layer that emerges as diffuse irradiance at the bottom of the layer, 
after scattering within it. 




£M 

Mo^o 


(16) 


Bihemispherical transmittance T 4 is the proportion of the incident diffuse irradiance that 
emerges from the bottom of the layer. 

1 

T* = 2^"/io(T, + r*)rf/Jo (17) 


L Non-Con— rvative Scattering 

When «< 1, directional-hemispherical transmittance is 

i-W> - 


# -V*0 


tills* 0 -S~*°) 




-fc(ciA«o+74) 


(18) 


2. Conner— lire Scattering 

When a = 1 directional-hemispherical transmittance is 

7* ( *“ l) (M 0 ) * 1 " UO +.^] = ^-x 

(l-e^( 7 ^o + 74) + 7lTo 


-3--1 

# -v*> 

,Mo 



(19) 
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